Post-Newtonian corrections to the motion of spinning bodies in NRGR 
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In this paper we include spin and multipole moment effects in the formalism used to describe the 
motion of extended objects recently introduced in 0. A suitable description for spinning bodies is 
developed and spin-orbit, spin-spin and quadrupole-spin Hamiltonians are found at leading order. 
The existence of tidal, as well as self induced finite size effects is shown, and the contribution to the 
Hamiltonian is calculated in the latter. It is shown that tidal deformations start formally at 0(v 6 ) 
and 0(v 10 ) for maximally rotating general and compact objects respectively, whereas self induced 
effects can show up at leading order. Agreement is found for the cases where the results are known. 
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I. INTRODUCTION 

In a recent paper an Effective Field Theory (EFT) 
of gravity for non spinning, spherically symmetric 
extended objects was introduced [y. Within the Post 
Newtonian (PN) framework this approach was coined 
NRGR (Non Relativistic General Relativity) due to 
its similarities with EFT approaches to non relativistic 
bound states in QED and QCD 0. However, the EFT 
formalism can be applied to a variety of scenarios, for 
instance the large small mass ratio case 0. NRGR 
is relevant for understanding the gravitational power 
spectra emitted by binary systems, an important class 
of candidate signals for gravitational wave observatories 
such as LIGO or VIRGO 0,E|. The formalism allows for 
a clean separation of the long wavelength gravitational 
dynamics from the details of the internal structure. 
This separation enables us to calculate corrections to all 
orders in the point particle approximation. Furthermore, 
it was shown that the ambiguities that plague 
the conventional PN calculations can be attributed to 
the presence of higher-dimensional worldline terms in 
the action whose coefficients encode the short distance 
structure of the particles. 

Building upon this idea, here we propose an extension 
of NRGR which allows for the inclusion of spin and mul- 
tipole moments. Spin in General Relativity (GR) has 
been considered previously in the literature from m any 
different points of view (see for instance 0, 0, 0, 
and references therein), and is argued to play an impor- 
tant role in binary inspiral, particularly for black holes 
[l3L H3| . Within the PN approximation spin effects for 
binary systems have been calculated using different tech- 
niques 

[30HHG3- Dealing with spinning objects 
in the point particle approximation inevitably entails 
running into divergent integrals as one does in the non 
spinning case. Regularization procedures, like Hadamard 
finite par t 1181 or considering contributions from differ- 
ent zones Il9| , were invoked when dealing with point- like 
sources |17| . However, it has been argued for instance for 
the proposal of |19|. that the formalism is "considerably 
complicated but it is inevitable that we have to adopt it to 
deal with divergences when we go to higher PN orders" 
[l6j| . As it has been repeatedly emphasized in that is 



not the case within an EFT approach. Here we will ex- 
plicitly see how a systematic, and consistent to all orders 
approach is translated into our case as well. 

The outline of the paper is as follows. In the first sec- 
tion we review NRGR for non-spinning spherically sym- 
metric objects highlighting the main results. Then, we 
will generalize the formalism to include internal angular 
as well as multipole moment degrees of freedom extend- 
ing the work of A. Hanson and T. Regge in the realm of 
special relativity [2(j ■ Afterwards we derive the power 
counting and Feynman rules of NRGR and calculate the 
leading spin-spin and spin-orbit potentials and show to 
reproduce known results [T(J, US 03 ■ A quadrupole-spin 
correction to the gravitational energy is obtained for the 
first time (to my knowledge) to leading order. The equiv- 
alence between different choices for the spin supplemen- 
tary condition is also shown. Finally we discuss the in- 
sertion of non-minimal terms in the worldline action and 
its relevance to renormalization. The existence of two 
types of finite size effects encapsulated in a new set of 
coefficients which can in principle be fixed by matching 
to the full theory, is predicted: those which have a renor- 
malization group (RG) flow and naturally represent tidal 
spin effects induced by the companion, and those that 
do not have a RG flow and represent self induced ef- 
fects, such as the spin induced quadrupole moment due 
to the proper rotation of the objects |2lL |22| . By power 
counting it is shown that companion induced tidal effects 
start formally at 3PN, and 5PN for maximally rotating 
general and compact objects respectively. Self induced 
effects can show up at leading order. Details are rel- 
egated to appendixes. We will study higher order PN 
corrections, the radiative energy loss, matching and new 
possible kinematic scenarios in future publications. 



II. NRGR 



In this section we will emphasize the main features of 
NRGR within the PN formalism, detailed calculation and 
further references can be found in the original proposal 
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A. Basic philosophy 

The traditional approach to the problem of motion 
was introduced by Fock [2j| who split it into two 
sub-problems. The internal problem, which consists of 
understanding the motion of each body around its center 
of mass, and the external problem which determines the 
motion of the centers of mass of each body. Decomposing 
the problem this way allows us to naturally separate 
scales and henceforth calculate in a more systematic 
fashion. The price one pays is the necessity of a matching 
procedure which relies either in comparing with the full 
theory, if known, or extracting unknown parameters 
from experiment. This method is now called "effective 
theory" , or EFT in the realm of Quantum Field Theory 
(QFT), and has been used to great success in many 
different branches of physics Q. While at first glance 
quantum field theoretical tools appear to introduce 
unnecessary machinery for classical calculations 1 , the 
power of the method will be shown to reside in two 
facts: It allows for the introduction of manifest power 
counting and naturally encapsulates divergences into 
text book renormalization procedures 2 . This means 
in addition that in the EFT it is straightforward to 
calculate the order at which a given term in the perturba- 
tive series first contributes to a given physical observable. 

Here we are going to tackle the problem of motion 
by treating gravity coupled to point particle sources 
as the classical limit of an EFT, i.e. the "tree level 
approximation" , within the PN formalism. Feynman 
diagrams will naturally show up as perturbative tech- 
niques to iteratively solve for the full Green functions of 
the theory. As it is known, GR coupled to distributional 
sources is not generically well defined due to its non 
linear character 3 [26j. This can be seen as a formal 
obstacle to the PN expansion for point particle sources 
in GR. Within an EFT paradigm this problem does 
not even arise since one is not claiming to construct 
a full description to be applicable to all regimes, but 
an effective theory which will mimic GR coupled to 
extended objects within its realm of applicability. In 
addition one can also argue that this EFT could be seen 
as the low energy regime of a quantum theory of gravity 
necessary to smear out point-like sources. 

The idea of describing low energy quantum gravity 
as an EFT is not new, for a review see j^]. What 
makes NRGR appealing is the uses of EFT to attack so 
called classical problems. QFT has proven to be useful 



1 QFT techniques have been recently used to calculate self-force 
effects in a curved spacetime background |24| . 

2 Also bear in mind that the classical solution is just the saddle 
point approximation to the path integral or what is known as 
the "tree level" approximation. 

3 I would like to thank Jorge Pullin for discussions on this point. 



with classical calculations, as in electromagnetic radia- 
tion where we can think of photons (QED) to calculate 
a power spectrum. Here we will use the same idea in- 
troducing "gravitons" as the quantum of the metric field 
which will allow us to calculate the gravitational poten- 
tial, from which the equations of motion (EOM) are de- 
rived, as well as gravitational radiation in a systematic 
fashion. 



B. Effective theory of extended objects 

The method of pj is based in the explicit separation of 
the relevant scales of the problem: the size of the objects 
r s (internal problem), the the size of the orbit r (exter- 
nal problem) and the natural radiation wavelength r/v, 
where v <C c is the relative velocity in the PN frame. 
Finite size effects are treated by the inclusion of a tower 
of new terms in the worldline action which are needed 
to regularize the theory 4 . For a non spinning spheri- 
cally symmetric particle the most general action consis- 
tent with the symmetries of GR is 



S= J {-m + c R R + c v v>*v v R IM/ + ...)dr, (1) 

where R^ v and R are the Ricci tensor and scalar 
respectively. The series, involving higher order Riemann 
type insertions, must be truncated within the desired 
accuracy to have any predictive power. The coefficients 
of each of these new terms can be determined by com- 
parison with the full theory. In this case the underlying 
theory is GR plus the internal equation of state of the 
objects. The beauty of this method is that, since these 
are 1-body properties, we can match using any relevant 
observable, for instance scattering processes, rather than 
solving the complete problem of motion explicitly. 

As it was shown in |lj the terms proportional to cr, cy 
are generated by logarithmic divergences of the point par- 
ticle approximation. However it is possible to show that 
they are unphysical in the sense that they can be removed 
from the effective action by field redefinition (f.r.) and no 
trace is left in observable quantities 0, ■ Nevertheless 
from here one concludes that not all divergences can be 
absorbed into the mass and new counterterms are nec- 
essary. Furthermore, at higher orders 5 it can be shown 
that (full Riemann dependent) finite size tidal effects are 
induced which can not be removed from the theory. We 
will see here that allowing the objects to spin also intro- 
duces new terms in the worldline action. For the sake of 



4 A similar approach can be found in l28l within the realm of 
tensor-scalar theories. However, renormalization issues as well 
as spin effects were left undiscussed. 

5 v 10 for non spinning particles. 
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completeness, and given that the same idea will be used 
here later on, we will sketch the reasoning. One starts by 
calculating the effective action pflj . 



k 
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(-k)T^ p) (k) + 
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where g^ — rj^ + , and /i M „ the graviton field. Let us 
concentrate on the contributions to the one point func- 
tion (k) . As it was shown in [lj within dimensional 
regularization techniques (dim. reg.), the logarithmic di- 
vergences in T£<.(k) can not be absorbed into the mass 
and a new counter-term of the form, 



T c 7(fc) = (2n)S(k-v) 



c R Wk 2 - Wk v ) + Jcv&W 



is therefore needed. It is straightforward to conclude 
from here the necessity of including two new terms in 
the effective action as shown in QJ. Within dim. reg. 
an arbitrary mass scale fj, associated to the substraction 
point at which the theory is renormalized is introduced. 
Given that the metric field does not pick any anomalous 
dimension at tree level we must have fj,dT[g^]/dfi = 0. 
Thus the explicit dependence on the subtraction scale 
fi must be cancelled by allowing the coefficients cr y to 
vary with scale. The theory therefore exhibits non-trivial 
classical RG scaling. As we are going to show here spin 
dependent finite size effects are predicted by similar ar- 
guments. 



C. The Post-Newtonian expansion 



TABLE I: NRGR power counting rules. 



The advantage of this redefinition is that now deriva- 
tives acting on any field in the EFT scale in the same 
way, <9 M ~ v/r, so it is easy to count powers of v coming 
from derivative interactions. 

The effective radiation NRGR Lagrangian, with the 
potential gravitons integrated out, can then be derived 
by computing the functional integral, 



(3) e-xp[iS NR GR[x a ,h]] = T>H^exp[iS[h+H,x a }+iS GF 



(8) 

where Sqf is a suitable gauge fixing term. Eq. (JSJ indi- 
cates that as far as the potential modes H^ u are con- 
cerned h^ v is just a slowly varying background field. 
To preserve gauge invariance of the effective action, we 
choose Sgf to be invariant under general coordinate 
transformations of the background metric g^ v (x) — rj^ v + 
hfj.„(x). This whole procedure is what is usually known 
as the "Background Field Method" , originally introduced 
by DeWitt 30] in canonical quantum gravity and used 
by t'Hooft and Veltman for the renormalization of gauge 
theories • By expanding the Einstein-Hilbert action 
using Q we can immediately read off Feynman rules . 
For potential gravitons, which we are going to represent 
by a dashed line, the propagator is given by, 



Once the internal scale is taken into account by the in- 
troduction of a series of new terms in the 1-body world- 
line action, the next scale we have to integrate out is 
the orbit scale. In order to do that we decompose the 
graviton field into two pieces, 



(4) 



where represents the off-shell potential gravitons, 
with 

1 v 
diH^ v ~ -H^v doH^ v ~ —H/iu, (5) 

and h^ v describes an on-shell radiation field 



(6) 



We can now further decompose by removing from 
it the large momentum fluctuations, 



H fiv (x) 



d 3 k 



_,ik-x 



(J) 



(H^x^H^piO)) = -(27r) 3 S 3 (k + q)—S(x )P^, a p, 

(9) 

where P^ v - a (3 = \ [VnaVvP + VvPVva ~ VnuVap]- Trie ra- 
diation gravitons, which will be represented by a curly 
line, have the usual spin 2 massless propagator. A wavy 
line will be used for the full propagator. We also need to 
consider mass insertions which will just provide a vertex 
interaction 



E 



-h 00 + h 0i v ai + i/ioov^ + 2 h ij v a 



where hoo,hoi,hij are evaluated on the 



(10) 

point parti- 



cle worldline (the leading order graviton-mass vertex is 
shown in fig. ^ . Following standard power counting pro- 
cedures we arrive to the scaling laws for the NRGR fields 
shown in table I 0, 0] . In the last column we have intro- 
duced mi = „„ V. the Planck mass and L = mvr the 
angular momentum. 

The effective action in (jSJ will be a function of the 
worldline particles (treated as external sources) and the 
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FIG. 1: Leading order mass vertex, the dashed line represents 
a potential graviton, whereas the wavy line stands for the full 
graviton propagator. 



with rju = (1, —1, —1, —1) the Minkowski metric and ej 
a set of / = 0..3 orthonormal basis vectors such that the 
tensor metric is diagonalized at the point. From now on 
capital Latin letters will denote internal indexes (notice 
that they transform in 50(3, 1) due the residual Lorentz 
invariance), the other conventions are as usual. Given a 
tetrad we can define its transport through the particle's 
worldline using Fermi- Walker ideas as |33j . 



radiation field which allows us to calculate the energy loss 
due to gravitational radiation as well as the gravitational 
binding potential from which the EOM are obtained. To 
get S e ff(x a ) we simply integrate out the radiation field 
h. S e ff(x a ) has a real part which represents the effective 
potential for the 2-body system 6 . It also has an imag- 
inary part that measures the total number of gravitons 
emitted by a given configuration {x 1 ^} over an arbitrarily 
large time T — > oo, 



1 If 

-lmS eff {x a ) = - dEdn 



d 2 T 
dEdtt'- 



(11) 



where dT is the differential rate for graviton emission 
from the binary system from which the power spectrum 
is computed. 



De'j 
~d\ 



(14) 



where V Q is the covariant derivative compatible with g, 
namely V a 9iiu = g^v.a — 0, and f7 MI/ is an antisymmetric 
tensor which therefore preserves l|12f) . One can invert the 
previous relation using (|12fl . 



f7fl Q 



u De aJ 
V e/3J- 



dr 



de aJ 
dr 



r^e^ju 7 ) r] IJ e 0I . 



(15) 



Notice that (|15|l implies the antisymmetry directly 
from g UKtl = 0. 



III. INTERNAL DEGREES OF FREEDOM I: 
SPINNING PARTICLE 

Here we will follow closely the ideas developed in [23] • 
We will start formulating a Lagrangian formalism to deal 
with internal angular as well as multipole moment de- 
grees of freedom which will enable us to describe a richer 
tensor structure. We will introduce the basic elements 
first, then we will construct the action and show how to 
reproduce Papapetrou equations for spinning particles in 
GR [2f|. The issue of constraints, the correct number 
of degrees of freedom and the angular- velocity/spin rela- 
tionship will be discussed at the end of the section. 



A. Basics 

Given a spacetime structure (g, M) we can always find 
at each point x €E M a coordinate system where the met- 
ric looks locally flat at the point. Such a transformation 
can be expressed as: 



V e I e J — 9 



(12) 
(13) 



The introduction of e^ 1 is equivalent to adding an ele- 
ment of 5*0(3, 1) to the worldline of the particle to de- 
scribe rotations [20| • Following these ideas we will there- 
fore construct an action in terms of the generalized coor- 
dinates and velocities (x^, u", e^, ejj. The number of de- 
grees of freedom in e^ 1 is 3 more than we need to describe 
3-rotations. We will see however that we can impose a 
set of kinematic constraints which will ensure the correct 
number. 



B. Action principle and EOM 

So far we have characterized the extra degrees of free- 
dom we need in order to construct a Lagrangian for the 
spinning particle. In the process to construct the ac- 
tion we will demand in addition to general covariance, 
internal Lorentz invariance as well as reparametrization 
invariance (RPI). This will naturally restrict ourselves 
to Lagrangians of the form L(x p ',u u ,il tJ ' u ). It is how- 
ever natural, instead of using W as coordinates to treat 
them as velocities of angular degrees of freedom which 
will lead us to a natural interpretation of spin. It is easy 
to see there are four different scalar quantities (neglecting 
parity violating terms) we can consider (schematically), 



Remember we are treating the worldline of the particles as exter- 
nal sources, namely x a = J, therefore S e ff is just the partition 
function e iS W/ =< 0|0 > J ~ e iC{jJ-)T ) witn r(j) the effective 
action for the sources as T — » oo |32| . 



di = U 

04 = £1575717 



(16) 
(17) 
(18) 
(19) 
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where contractions are made with the space-time 
metric g^. Using these quantities our Lagrangian 
will be in principle a general expression of the form 
L(ai, ci2, (I3, 04). We will neglect multipole moments 
throughout this section the inclusion of which will be 
studied later on. The objects will be therefore considered 
symmetric with respect to their rotational axis. 

In order to introduce the idea of spin we will define the 
antisymmetric tensor S^ u and momentum p* 1 by, 



ativistic limit 20]. From these definitions we will have, 



P 



= -2u c 



dL 
dai 



up ' 



,dL 
da 3 



(21) 



f)T r)T 

-40^- 2(u^ A u A - ^uxItt- 

oa 2 oa 3 



00,4 



dL 



(22) 



1 



sl = -p^Su^ - -s^sn^, (20) 

where the minus sign corresponds to the correct non rel- 



The variation of the action consists of two pieces. Let 
us concentrate first in the tetrad part. Using the defini- 
tion of spin we will have to deal with, 



J 



SS 



dTS al3 dn a fj = - I dr(- 



DS af3 

~l>t~ 



-&Kc 



De Ka 
Dt 



S 



Dt 



6eZ. 



r 



(23) 



The equation of motion can be directly read from the 
above expression, multiplying by e Ktl we get (using l |15t ). 



DS^ 
Dt 



p»u v 



(24) 



where the last equality follows from I|21I22|I 20]. Notice 
that we have not specified a Lagrangian up to this stage. 
It is easy to see from (|24ll that, 



since it is covariant. We therefore recognize in (|24I26|) 
the well known Papapetrou equations |25| . Remarkably, 
although in terms of the tetrad these equations depend 
on the choice of Lagrangian, as a function of spin and mo- 
mentum the evolution equations are action-independent 
as far as curvature terms are not inserted. Including cur- 
vature terms in the effective action will be relevant to 
introduce finite size effects and will modify these equa- 
tions for extended objects. 



DS 



1. j 



D(S 



Dt 



Dt 



= 0, 



C. Constraints and angular- velocity /spin 
(25) relationship 



which shows that spin projected with respect to the 
frame remains constant. In addition one can also show 
that the scalar S 2 = \S^ P ' is conserved. As a further 
property it is also instructive to notice that 5 M " S*^ v 

if 



is also a constant of the motion, where S 



In order the get the 8x piece of 8S, a shortcut can be 
taken by going to a locally flat coordinate system where 
the connection terms are zero at the point. Promoting 
the derivatives to covariant ones in this frame we will end 
up with, 



Dt 



-S P (Pa/3cr,7 ^a/3j,<r)u 



(26) 



where we have used the form of the Riemann tensor in 
a locally flat coordinate system. Written this way we 
can promote now the equation to all reference systems 



In order to describe the correct number of degrees 
of freedom we need to add a set of constraints to the 
EOM I|24I26|) . A well defined angular- velocity/spin rela- 
tionship is also necessary to extend our power counting 
rules to the spinning case. We will show here that both 
features are related. Here we will closely follow [2(|, to 
which we refer the reader for details, other approaches 
may be found in [34|. This section relies on a basic 
knowledge of constrained systems, for further details see 

MM. 



It is natural to impose the following (covariant) con- 
straints in the space of solutions, 



U" = S^p v sa 0, 



(27) 



where just three of the four components are indepen- 
dent, and stands for weakly vanishing 40]. This 
set of constraints are second class, namely they have 
non vanishing Poisson bracket among themselves, and 
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therefore reduce the number of degrees of freedom 
from 6 to 3 SO(3) parameters as expected. It can be 
shown in addition there is a Lagrangian from which Ij27(l 
kinematically follows pcj . 

We need to guarantee the constraints in Ij27(l are pre- 
served upon evolution. It is however possible to show 
from (|24l2tj|) that ^f- w will be satisfied provided 

p a = mu a -^-Rp upa S al3 S pa u v , (28) 

with m 2 {S 2 ) = p 2 defined by I|21I22|I . This also means 
that the difference between p v jm and u u is higher order 
in the PN expansion and we can consider S^ v u v = as 
well as ds , — to leading order. 

It is possible to show that |2ZJ) implies C\ = S<*"S* U w 
0, which is a first class constraint. There is therefore, in 
addition to RPI (C*2 = p 2 — m(S 2 ) « 0), a gauge freedom 
which can be attributed to the choice of the temporal 
vector of the tetrad e°, and a sensible choice of gauge is 
then = Cq — p^/m w [2(j. This gauge, jointly with 
(|27|l . also translates into a choice of center of mass of the 
object |ToL Il4|. and implies as well VL tiV p v = from 
which we get 

~ Snv — RnuapS ^ + RRS S S + (29) 

2m 

where we have used I|27l28fl . We can indeed obtain the 
angular- velocity /spin relationship by matching the evo- 
lution equation for the tetrad in a Minkowski background 
obtaining (see appendix A), 

S " V = (^172^ + h R ^ a " + (30) 

with I the moment of inertia, and a particular La- 
grangian is chosen to ensure l|27|l pp| . The main results 
of this section are therefore equations (|28I3U[) . from which 
we conclude that in a theory where (|27|l is kinematically 
imposed spin and angular-velocity are naturally related 
and proportional in flat space. Within an EFT approach, 
these relationships are all we need to construct the NRGR 
extension for spinning bodies. 

IV. NRGR FOR SPINNING BODIES 

A. Power counting 

The power counting rules in NRGR have been devel- 
oped in p] . Here we are going to extend them to include 
spin degrees of freedom. As we shall show the only nec- 
essary change is the inclusion of spin insertions at the 
vertices. 



First of all notice that, from the constraints at leading 
order, 

S» v u v = -» S 3 ° = S lk u k , (31) 

which implies a suppression of the temporal components 
with respect to the spatial ones. In other words spin 
is represented by a 3- vector, S — he in the rest 

frame of the particle as expected. We will concentrate 
here in compact objects like neutron stars or black holes 
where the natural length scale can be taken to be their 
Schwarzschild radius r s ~ Gm (for general objects see 
appendix B), and hence a momentum of inertia scaling 
'as I ^ m 3 /nip. For the spin angular momentum we will 
have, 

S = ltd = 1—^- ~ mv ro tr s < mr s ~ Lv. (32) 

We therefore see that spin gets suppressed with respect 
to the orbital angular momentum, even for the maximally 
rotating case (v rot = 1). We can also assume a different 
kinematical configuration with the particles co-rotating, 

namely = -. In the former S ~ m " rs = ~ Lv 4 . 
We will generally power count spin as 

S - Lv s , (33) 

with s = 1 and s = 4 in the maximally rotating and 
co-rotating scenario respectively. 

To obtain the scaling laws we have assumed the usual 
proportionality at leading order between spin and an- 
gular velocity which was obtained from iJ5U|) . As usual 
subleading scalings will be naturally taken into account 
by the insertion of higher order terms in the worldlinc 
action. What we have learnt here is that spin effects are 
in any case subleading in the PN expansion, and the scal- 
ing laws developed in still hold and spin contributions 
can be treated as a perturbation. 

B. Feynman rules: spin-graviton vertex 

To construct the effective theory for gravitons we need 
to expand the metric around a Minkowski background, 
namely g a p = r\ a p + . There is however a subtle point 
in doing this given that l|13fl leads to, 

fluent = V„u + ^, (34) 

where A Ia is an element of the Lorentz group. We will 
also need Se^ 1 which is defined through the inverse metric 
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FIG. 2: Leading order spin-graviton vertex interaction, 
blow represents an spin insertion. 



The 



will differ in the spin/angular- velocity relationships and 
might not lead kinematically to (|27|) . The physics will be 
invariant once these differences are taken into account. 

To calculate in the EFT we need to match (|38I39[I into 
NRGR using the power counting rules developed in p| 
plus the spin insertions. Up to 2PN, for maximally ro- 
tating compact objects we will get for the spinning part 
of the NRGR Lagrangian, 



h?" = ri^'v vv '(-hn>v>+h%h au ' -...), (35) 

where indexes are raised with rf^ . We will therefore 
have, 



2m, 



-rf v h va K al + 



(36) 



L 



NRGR 


1 


1PN — 


2m p 


NRGR 


1 


'1.5PN — 


2m p 


NRGR. 


1 


'2PN — 


2m p 




1 


+ 


Amp 1 



i0,0<- 

k: 



(40) 



[Hij, fc S* V + H 0Qtk S° k ) , (41) 
(H 0i . k S ok ui + H l0 . S M ) 



_V<(//f// ::A ., //j //,,,). (42) 



One can immediately see how to proceed by compari- 
son with what it has been done in |20| within flat space 
where the angular velocity was defined as, 



dr 



(37) 



with A Ml/ describing the rotation of the particle, and 
M stands for Minkowski. Expanding the action using 
(j34l36[l we will thus obtain the spin-graviton interaction 
to all orders in a flat space background as a function of 
the graviton field and (|37|) (see appendix C for details). 
To leading order in the weak field expansion(see fig. 0), 



u = 



2m 



■a/3 



~h a ~/ fiS 



(38) 



The procedure follows systematically as shown in ap- 
pendix C. 

It is an useful exercise to check the gauge invariance of 
l|38jl. or in other words to obtain the leading stress energy 
tensor. It is straightforward to calculate = ^t^- a t 

sv = V g ettm S> 



Tfi = -~dp (S^u u + fi^u") , (43) 



which agrees at zero order with the original proposal of 
Dixon |3(| and also Bailey and Israel (BI) [33 (see ap- 
pendix D). 

The Ward identity, 



where S 1 ^ = — at g^ v = r)^ (we will drop the 
M from now on). The expression in l|38[) is remarkably 
action independent if written in terms of spin and the 
graviton field. The Lagrangian dependence enters in the 
so far unknown function S(Q). This conclusion however 
just applies to the leading term and a choice of action is 
necessary to obtain the Feynman rules to all orders. At 
next to leading order in the weak gravity limit we will 
have, 



Li 



4m£ 



(39) 



where a particular action has been chosen to ensure (|27H . 
A different choice of Lagrangian will imply different Feyn- 
man rules. However, bear in mind that different actions 



(44) 



is therefore obeyed since spin is constant to leading order 
in the PN expansion. 



C. Leading order graviton exchange 

Our goal from now on is to calculate the leading or- 
der piece (one graviton exchange) of the potential energy 
due to spin-orbit (fig. |3J) and spin-spin (fig. 0} couplings 
coming from 140I41|) 7 . The leading order spin-orbit con- 
tribution is the sum of two pieces, 
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7 Self energy terms are not considered since they yield scaleless 
integrals. 



J 



fig.®-- 



— im 2 

Am 2 



dtdt' 



d 3 p 
(2^f 



-ip(x(t)-y(t')) 



P 2 



-S(t - t')P 0e - ai I S( a uJ(t)u 2 (2 - 5 e Q ) 



(45) 



We need to distinguish two different cases: the tempo- 
ral and spatial derivative. The temporal derivative will 
just hit 5(t — t') and can be integrated by part bringing 
down a velocity factor. To leading order we would need 
to consider e = 7 = 0. However, Poo;aO = unless a = 0, 
which leads to a term proportional to S 00 — 0. There 
is then no contribution from the temporal derivative and 
we just need to concentrate in the spatial part and the 
terms, 



-im 2 G N dj f dt T ——^-——P 0e . ai S{ a 
J !*(*) -»(*)! 



u\{t)ul{t){2-5l). 

(46) 

Notice that we have three possible contributions, one 
coming from e = 0, a — i, 7 = k, another where e = 
I, a = k,^f = and finally e = a = 7 = 0. The latter 
looks at first as a v° piece, however this is misleading 
since our spin choice implies S^ — S^u 1 . Adding all the 
terms one gets, 



v" 
b) 



FIG. 3: Leading order spin-orbit interaction. Diagram a) 
takes into account a v 1 mass insertion. Diagram b) corre- 
spond to the v 3 spin-graviton vertex. 



fig®= i / dt 



\m-m\- 



^(n x ui) ■ S\ — (n x 112) ■ Sij 



(47) 

where n is the unit vector in the (x — y) direction and 
we have used — e k: > l S l . Joining the mirror image we 
will end up with, 



V so 



2Gn 



/i(n x v) 



1 + ^)S 2 

m 2 



mi 



for the spin-orbit potential, where /i is the reduced mass, 
r = \x(t) — y(t)\ and v = u\ — v.2- 

Let us now consider the spin-spin interaction. The 
leading order contribution is (see fig-EJ, 



' B f) I rlf r ip(S(t)-v(t)) 1 p gjkgfk' 



(27T) 3 



V 2 



(49) 



Using Pj'o-jo — h is straightforward to show 



V ss 



G 



N 



Si ■ S 2 — 3Si ■ nS 2 ■ ft 



for the spin-spin binding potential. It is easy to see by 
power counting that Vsodt ~ Lv 3 and Vssdt ~ Lv 4 , ef- 
fectively 1.5PN and 2PN for maximally rotating compact 
objects. 



FIG. 4: Leading order spin-spin interaction. 



By comparison with the results in |lf)j | , it is immediate 
to notice there is a mismatch in the spin-orbit contribu- 
tion H48fl , which can be traced back to the choice of spin 
sup p lementar y condition in 127|) . As it was noticed in 
[nj,[l2,[Ll,[l]| this discrepancy is associated to the choice 
of center of mass of each body. It can be shown there is a 
coordinate transformation that relates the center of mass 
choice which follows from l|27|l and the so called baryonic 
coordinates (implicitly used in jlflj). where the center of 
mass is defined through the baryonic density, and one 
has S M = iS^ui. Had we calculated the spin-orbit term 
within the baryonic condition we would obtain, 



(50) vso 



2G 



N 



3mi 



377J2 



(51) 

in complete agreement with the result in [lOj. The lead- 
ing order spin-spin interaction does not get affected by 
this new choice. 
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D. EOM 

Even though 1|51|) is in total agreement with the spin- 
orbit potential in baryonic coordinates, the calculation 
in the covariant approach does not reproduce the well 
known fact that the generalized Lagrangian from which 
the EOM are derived turns out to be acceleration de- 
pendent 0, E IH Hil. I ndeed. igHjl reproduces the 
gravitational potential in [13j up to this acceleration de- 
pendent piece which does not follow from a graviton ex- 
change. As we shall show in Appendix E, the solution 
to this puzzle lies on the fact that a non canonical alge- 



bra develops which naturally reconciles both approaches. 
Instead of following that path here it is instructive to re- 
mark there is a coordinate transformation which leads to 
a canonical structure. Not surprisingly this map trans- 
forms the covariant choice into the baryonic one, where 
it has been explicitly shown there is no need for an accel- 
eration dependent piece in the action. We can therefore 
proceed from the potentials in l|5UI51|l and the standard 
Euler-Lagrange formalism to obtain the EOM within the 
baryonic supplementary condition. It can be easily shown 
that they are given by (in relative coordinates), 



a = <2l — 0,2 

aso 

ass 

dSi 
~dt 



GnM _ _ 
2 — n + aso + ass 



G 



N 



(3n(n x v) ■ \ + 2v x x + 3ri ■ v(n x x)) 



3G 



N 



ft(Si ■ S 2 - 5S 2 ■ nS 2 ■ ft) + S 2 (n • Si) + Si(n ■ S 2 ) 



^-)+S 1 xS 2 + 'S{n-S 2 )nxS 1 J; -^ = 1^2, 



(52) 
(53) 
(54) 

(55) 



r 



where we have introduced M — mi + m 2l x = 
( 2 + ill) Si + (2 + S 2 and L = »rft x v. 

From the symmetries of the action we can directly con- 
struct the conserved quantities, in particular the energy. 
First of all notice that the spin-orbit force does not do 
any work, namely aso • v — 0, from which we conclude 
that the conserved energy is nothing but 



E = -tf 



G N M/i 



V ss - 



(56) 



In order to compare with the results in [l3| within the 
covariant supplementary condition, we restrict ourselves 
now to the case of nearly circular orbits to express (|56|l in 
terms of the orbital angular frequency u> and spin. Taking 
an angular average for all quantities we obtain from 15211 , 



M 1 ' 3 

LO 2 / 3 



' X - 77 ; rR /3 I Si ■ S 2 -3l ■ Sil ■ S 2 



3M 



2 11M 5 / 3 

uj 4 ' 3 



E(u,S) = --{Muj) 2 ' 3 ( 1 + T-x+ - . 



S\ • S2 — 3/ • S\l • S e £ 



(57) 
(58) 



with I the unit vector in the L-direction. We therefore 
conclude that energy as a function of spin and angular 
frequency in l|58|) matches that of |13l Il5j independently 
of the spin supplementary condition. As a final comment 
let us remark that our results also agree with those of 
Buonanno et al. that appeared after we had com- 
pleted our work, where a similar procedure is advocated 
within baryonic coordinates. However, the covariant spin 
supplementary condition and finite size effects are left 



undiscussed in [38| . 



V. INTERNAL DEGREES OF FREEDOM II: 
PERMANENT MULTIPOLE MOMENTS 

By now it should be easy to visualize how are we going 
to include multipole moments in terms of the fields. 
Let us assume for instance the particle has an intrinsic, 
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FIG. 5: Leading order quadrupole-spin one graviton ex- 
change. The black square represents a quadrupole insertion. 

permanent quadrupole moment Q IJ . In order to couple 
it to the gravitational field the following term can be 
introduced, 

J dTRp a fre$ePjQ IJ <u a ui '. (59) 

Notice in fact this is just a generalization of what it 
has been done in [l|. In fact, it can be shown (see ap- 



for the quadrupole-spin potential within the covariant 
spin condition, where we have used Q y = Qi % and the 
Euclidian metric (<5 y ) to raise and lower indexes. A 
similar expression is obtained from the mirror image 
1 <-» 2. It is easy to show it corresponds to a 3.5PN con- 
tribution for maximally rotating neutron stars or black 
holes. Notice that for a spherically symmetric object 
(Q y ~ <5 y ) Eqs — > as one would have guessed 8 . There- 
fore, the coupling is effectively to the traceless piece of the 
quadrupole. Higher order multipole moments are easily 
handled by similar procedures. 

A. Quadrupole radiation 

It is instructive to notice that (|59|l will directly lead 
to the well known quadrupole radiation formula. The 
leading order piece will be of the form (for on-shell gravi- 
tons), 

-^-RiOjoQTF' (61) 

where TF stands for the traceless piece. By calculating 
the imaginary part of fig. we can immediately obtain 
(we skip details which can be found for an identical cal- 
culation in 0|), 



This is reminiscent of the vanishing of the cy contributions in 

u 



pendix C) this quadrupole term is naturally obtained if 
we include the non-spherical contribution of the tensor 
of inertia in the spin part of the Lagrangian. 

From the gauge fixing condition eg ss p M /m we imme- 
diately see that the Q 01 components do not contribute at 
all given that replacing in i|59[) give rise to vanishing 
terms. We can therefore, as expected, concentrate just 
in the spatial components. 

It is easy to see that (|59|) will naturally reproduce the 

quadrupole gravitational energy piece, < ^' jX s x ; in the po- 
tential. In order to obtain a non trivial contribution 
we calculate a correction to the binding energy due to 
a quadrupole-spin interaction to leading order. After 
matching H59|) into NRGR and using the expression for 
the Riemann tensor in the weak gravity approximation, 
the one potential graviton exchange in fig. [3] gives, 



(60) 

I 



1 f f{i\r 

,m/ '* H -^/(Sp2M k4|<? -< |k|)|2 ' (62 » 

from which the power radiated follows, 

S 1 p DO 

= Q'ri Q'ri ■ (63) 

whit dots as time derivatives and the bracket represent- 
ing time averaging. This is the celebrated quadrupole 
radiation formula. 

VI. DIVERGENCES, NON MINIMAL 
INSERTIONS AND FINITE SIZE EFFECTS FOR 
SPINNING BODIES 

We are going to discuss here the appearance of diver- 
gences and their consequent renormalization. This will 
lead us to the study of higher order terms in the worldlinc 
action and their Wilson coefficients, which will encode 
the information about the internal structure of the body. 
In addition to terms coming from logarithmic UV diver- 
gences, we will encounter power law divergences whose 
associated Wilson coefficients do not have scale depen- 
dence. This distinction will turn out to be connected 
with tidal deformations vs. self induced effects as well 
shall see. 



V QS = 



3G 



iV 



(Qil - 5Q\ k n t n k )n ■ (u 2 x S 2 ) + Q{ k n k (j2u 2 + uj.) x S 2 J + Q l { (u Xj - 5(ui • n) nj ) (n x S 2 ) t 
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FIG. 6: Leading order diagram whose imaginary part give ? 
rise to the quadrupole radiation power spectrum. The black <T 
boxes are quadrupole insertions and the curly propagator a <T 
radiation graviton. < - 

FIG. 8: Three graviton contribution to the one point function 
^ in the effective action with two spin insertions. 




FIG. 7: Two graviton contribution to the one point function 
in effective action with a single spin insertion. 



B. Non minimal insertions I: Self Induced effects 

We will consider here terms which are not total deriva- 
tives and can not be removed by f.r. Let us proceed sys- 
tematically. Let us start with terms linear in Riemann 
and no further derivatives acting on external fields. The 
first non zero terms we can construct are, 



A. A cursory first look 

Let us study the one point function in the effective ac- 
tion with spin insertions. Let us start with the diagram 
shown in fig. The spin-graviton Feynman rules derived 
from (|38|l differs from mass insertions in two main points: 
its tensor structure and its dependence on the graviton 
momentum. Fig. [7| then contributes for potential gravi- 
tons terms proportional to, 



qk,q 2 ,k 2 



(2tt)3 Vq 2 (q + k) 2 



(64) 



where k is the external graviton momentum. None of 
these integrals are logarithmically divergent and there- 
fore can be absorbed as pure counterterms in the original 
Lagrangian. 

Let us concentrate now on the divergent piece coming 
from diagrams like in fig. [SJ It can be shown that this 
diagram contain terms such as (in d-dimensions) , 



7(k) 



p rf'-iq (q-k)(p-k) 

{2n) d - 1 {2n) d - 1 q 2 p 2 (q + p + k) 



(65) 



as well as integrals with q • p in the numerator. These 
integrals contain power as well as logarithmic UV diver- 
gences. It is clear that these divergences can not be ab- 
sorbed into the original Lagrangian since they involve in 
principle higher order derivatives of the metric. By gen- 
eral covariance and parity conservation, there is a limited 
set of non zero terms one can build up with the right 
structure to cancel the previous divergences. In what 
follows we will study a possible set of new insertions in 
the worldline action which are generated by renormaliza- 
tion. 



Ors 2 



&RS 2 



Ors 2 



m p 

^-R a 8uuS a ^S^u". 



(66) 



It is possible to show that they are physically equiva- 
lent, namely, they are proportional up to f.r. removable 
terms. By simple inspection it is easy to see they are also 
similar to the quadrupole moment insertion in (|59|I .It is 
therefore natural to associate these operators with self 
induced quadrupole effects rather than tidal deforma- 
tions. It is straightforward to show, after matching into 
NRGR for potential gravitons, that O rs2 



y/Lv 



2s+2 



with C Rg2 ~ 1/m. This immediately tell us O rs2 can 
not be generated from a logarithmic UV divergence. It 
is possible to show nonetheless that these terms can be 
generated from power law divergences (see appendix G). 
At 2(5)PN these terms generate a gravitational poten- 
tial for a maximally rotating (co-rotating) neutron star 
or black hole A coupled to a non spinning one B of the 
form, 



v f 



s 2 o 



= C 



RS 2 (A) 2r3 



with C%p {A) 



(• 



3(S~ A ■ n) 2 - S A ■ S A ) , (67) 



1^2 



\ C rs 2 )a- a spin- 
ning particle will tend to deform and therefore generate 
multipole (mass) moments which will thus produce 
a binding energy term e quiv alent, as in this case, to 
a quadrupole interaction 9 |2jJ. It is well known that 



9 It will in principle vary in time and henceforth radiate. This 
effect will be naturally taken into account similarly as we did for 
the quadrupole. 
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rotating black holes, or neutron stars, have a quadrupole 



moment given by Qbh = —aS /m (G 



1), with 



DpDeR 



S°»S pa u v v?\ D a D 2 R 0pp , e S pa S'S^u' J ' 



m, S the mass and spin respectively. For a black hole 
a = 1 |39j]. for neutron stars a ranges between 4 and 
8 depending on the equation of state of the neutron 
star matter |22^ . This will give us an straightforward 
matching for Cjjg2^<n in GSt which is consistent by 
dimensional analysis with what we expect from nat- 
uralness arguments. Furthermore, these coefficients 
contribute to the one point function, and thus will 
show up in the metric solution for a rotating neutron 
star or black hole. For the case of a black hole, the 
Kerr-Newman spacetime does not have any logarithmic 
dependence and therefore, every coefficient associated 
to a non f.r. removable term 10 , which contributes to 
the one point function must be scale independent. 
Similarly for neutron stars. This provides a natural 
characterization for self induced effects. Tidal effects in 
the other hand will be then associated to coefficients 
which do not contribute to the one point function, and 
moreover are scale dependent. Tidally induced effects 
will be therefore naturally generated by logarithmic UV 
divergences. 

Once self induced spin-multipole moments are included 
it is no longer necessary to introduce a non dynamical 
permanent multipole. Adding more spin insertions with- 
out derivatives will have the same type of behavior we 
encountered above, namely the Wilson coefficients will 
scale with negative powers of the mass and therefore they 
can not be generated from logarithmic UV divergences. 
It is indeed possible to show that we can in principle 
hook together n spin tensors leading to terms scaling as 
^/Zyi+nis-i) after matching into NRGR. Given that n 
could be any number of spin insertions, it appears as if 
we will have no predictive power for the case s = 1 . How- 
ever, it can be shown that terms with a large number of 
spin insertions can be rewritten in terms of interactions 
with no more than four spin insertions [2(j. The case of 
three and four spin insertions does not modify the leading 
order expression in (|67|) . 



C. Non minimal insertions II: Tidal deformations 



We will concentrate here in tidal effects and therefore 
in those terms generated by logarithmic UV divergences. 
To lowest order it can be shown we will have D 2 O l RS2 
(with i = 1,2,3) coming from diagrams like fig. HI We 
will generically denote their Wilson coefficients by Cjyi . 
Given that these expressions posses different tensorial 
structure the RG flow will naturally decouple. By di- 
mensional analysis it is easy to conclude that 



dC r 



df i 



m 



(69) 



As it was pointed out before, these new insertions will 
not contribute to the one point function 11 . However, 
they will in principle be observable for more complicated 
ambient metric, such as the field produced by a binary 
companion. 

Let us power count this effect. After matching into 
NRGR for potential gravitons we will get, to leading or- 
der, 



6+2s 



(70) 



which would make it a 4PN contribution for maximally 
rotating compact objects. A careful inspection shows 
however that the leading piece from these terms goes as 
derivatives of 8(x\ — x-z), which is a contact interaction 
(the k 2 piece cancels the propagator). As a consequence, 
the first long range interaction coming from Cjp., and 
therefore the lowest companion induced tidal effect, will 
scale as Lv 8+2s , a 5PN contribution for maximally ro- 
tating compact objects (formally at 3PN as shown in 
appendix B). At this order new terms will also start to 
contribute (for instance the third and fourth expressions 
in (|68"|) . See appendix G for details). The reasoning in 
previous sections can be easily extended to the n-point 
function and higher Riemann insertions. 



Other type of worldline insertions we could in principle 
generate are those having derivatives of the Riemann ten- 
sor and more spin insertions. We will need to introduce 
terms like, 

D e R a ^S^S uf3 u a ; B 2 R llvaSi S^S a a v v u li - (68) 



As it was noticed in Q, f.r. removable terms can in principle 
show up in the one point function. However, they can be washed 
away by a coordinate transformation. Here we will concentrate 
in terms which are not f.r. removable. 



VII. CONCLUSIONS 

In this paper we have extended the formalism initially 
proposed in pj to include internal degrees of freedom 
like spin as well as multipole moments. As a first step 
we have developed in a suitable fashion the description 
of spinning bodies in GR to include a richer tensor 
structure extending the previous work done in the realm 
of special relativity by A. Hanson and T. Regge |20|. We 



11 For instance D 2 O i „ r , = on shell. 

H.S A 
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have shown that a self consistent action principle can be 
implemented and Papapetrou equations |25| recovered. 
Permanent multipole moments are naturally introduced 
by adding new degrees of freedom in the worldline action. 
Using this formalism we have extended NRGR, its power 
counting and Feynman rules with which we have repro- 
duced the well known spin-spin and spin-orbit effects at 
leading order [13J. A quadrupole-spin correction to the 
binding energy was obtained for the first time (to my 
knowledge) as well as the quadrupole radiation formula 
recovered. The equivalence between different choices 
for the spin supplementary condition was explicitly 
shown. We have shown afterwards the appearance of 
divergences at higher orders in the PN expansion and 
its consequent regularization. The type of divergences 
are twofold: logarithmic and power law UV divergences. 
This distinction was shown to be associated to tidal 
vs. self induced effects. Rcnormalization through the 
insertion of non minimal terms in the effective action 
was implemented and the RG flow obtained. A finite size 
cutoff was invoked in the case of power law divergences 
and its respective Wilson coefficients set by naturalness. 
In the EFT spirit it is likely that all terms which are 
consistent with the symmetries will contribute to the 
effective action. In fact, self induced spin effects are 
naturally expected and the lack of scale dependence just 
responds to the fact that it is a 1-body effect on itself due 
to its proper rotation which does not get renormalized 12 . 
A partial matching into the full theory was accomplished 
by comparison with known results |2ll |22( . Self induced 
effects could in principle appear at leading orders in the 
PN expansion in the case of maximally rotating compact 
bodies, for which tidal deformations were shown to first 
appear at 5PN, although formally at 3PN for general 
objects. Within the power of the EFT, most of the 
conclusions are based on dimensional grounds without 
detailed calculations. 



Several aspects remain still to be worked out. In ad- 
dition to the matching calculation and the issue of fi- 
nite size effects, higher order corrections are yet to be 
obtained even though the formalism is already set and 
just computational work is needed. Including dynamical 
properties for multipole moments as well as back reaction 
effects is also to be worked out. Moreover, new kinemat- 
ical scenarios, like a 3-body system and the large small 
mass ratio case, are currently under study. All these is- 
sues, including the radiative energy loss due to spin will 
be covered in forthcoming publications. 



12 Formally speaking, self induced effects do not get renormalized as 
a consequence of the fact that they are derived from the coupling 
to the conserved stress energy and the metric field does not get 
renormalized classically. 
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APPENDIX A: ANGULAR- VELOCITY/SPIN 
RELATIONSHIP 

By RPI we know the theory has vanishing Hamiltonian 
and dynamics is generated by the constraints Ci,C2. It 
can be also shown that the Lagrange multiplier associ- 
ated to Ci has been set to zero by the condition ip ss 
[2£j ■ Using the Hamiltonian equations for the tetrad and 
position we have in the realm of special relativity, 



2m 



[el^C 2 ] pb ^2U'(S 2 )S^ei 



(Al) 



where £ is a Lagrange multiplier, [, ] p t, stands for the Pois- 
son bracket and f(S 2 ) = m 2 (S 2 ). By comparison with 
(|29[1 we conclude, by matching to the zero curvature case, 



(ua) i /a /W (5 



1 

2m 



R^ a pS af3 + ...)■ (A2) 



The Lagrangian dependence of this expression is en- 
coded in the function f(S 2 ) defined bv_ (|2H22l) . It is 
possible now to construct a Lagrangian (L) using all the 
freedom we showed previously, that will ensure Ij27|l kinc- 
matically j2lj. Such a procedure is therefore preferred 
given that the unphysical degrees of freedom are cut off 
kinematically rather than cut by hand. As it has been 
shown in 20] L is however not unique. There is still 
a remnant freedom of the form f'(S 2 ) ~ A, with A a 
constant 13 . We can henceforth set A in order to recover 
the well known relationship between angular- velocity and 
spin in flat space, namely S ~ with / the moment of 
inertia. One then solves for S 1 *" in (|A2|) order by order 

to get dsnj ■ 

One could still argue that the angular- velocity/spin 
relation should be obtained directly from Sm — (§n) 



Regge trajectories are of constant slope. 
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rather than using the EOM as we did. One should how- 
ever bear in mind that dynamics naturally help us to 
power count within an EFT approach. Higher order cor- 
rections are taken into account by insertions in the world- 
line action 0. 

APPENDIX B: FORMAL POWER COUNTING 

Here we will comment on the power counting from a 
formal point of view without assuming any specific prop- 
erties of the objects. This will introduce new parameters 
which should be adjusted depending on the constituents. 
For an object of characteristic length R the spin magni- 
tude for a rotating velocity v rot is S ~ L^'^ L 1 formally 
a non perturbative effects. For maximally rotating, and 
co-rotating, bodies one gets S ~ L^- and S ~ L4 re- 
spectively. Nevertheless, it is naturally expected that 
e = R/r <C 1 and a new perturbative parameter is intro- 
duced (e ~ v 2 for neutron stars or black holes) . It is easy 
to show that the leading spin-orbit effect for maximally 
rotating objects scales as veL, a subleading contribution 
formally at 0.5PN. The first long range contribution from 
Cjj2 to the effective action will now scale as Lv e e 2 for 
Vrot = 1. A 3PN effect, effectively at 4PN for e ~ v. 
This brings hope to potentially observe these effects in 
the future. 



APPENDIX C: GOING TO ALL ORDERS 

As we pointed out the unphysical states can be washed 
away kinematically by l|27|l if a suitable Lagrangian (L) 
is chosen |2£|. Moreover, the leading order spin-graviton 
vertex was shown to be Lagrangian independent to 
leading order. This however can no be translated to 
higher orders. Here we will show how to proceed to 
obtain the Feynman rules to all orders. 

By local Lorentz invariance and general covariance 
we know L is a function of the metric and angular- 
velocity. We can rewrite L = L(£l IJ ,i] IJ ) which shrinks 
to Vt IJ = ejCjfl^ all the metric dependence. Within an 
EFT framework the explicit form of the Lagrangian given 
in |2fl | is not necessary, since we can always obtain its 
NRGR counterpart by expanding L around a Minkowski 
background, 

(Cl) 



where 5n IJ (5e, h) = n IJ ~fl^, and Sl$ defined by 
Using that Sm = I^m on shell, the NRGR Lagrangian 
turns out to be (schematically), 



l = l(q m ) - -s M sn - -snsn. (C2) 

By expanding SCl(Se, h) in (|C2(1 we will therefore 
generate the spin-graviton vertices to all orders in the 
weak field limit 14 . The next step to construct the EFT 
is to match into NRGR using the power counting rules 
thus far developed y|. That is an straightforward task. 
The terms in the NRGR Lagrangian for maximally 
rotating compact objects are shown in 140I41I42|) up to 
2PN. 

Notice also that the second piece in (IC2|) generates 
contributions which are not explicitly spin dependent, 
although the coupling is proportional to the moment of 
inertia. For spherically symmetric objects we will have 
for instance a term of the form, 



\ Y^T^u a uP, (C3) 

which can be easily shown to be proportional to R^u^u" 
and henceforth f.r. removable. For non-spherical bodies 
we will get, 

R^ a pe K ejI KJ u»u a , (C4) 

with I KJ the inertia tensor defined by I KJ = 
J2 p m p S KJ — XpXp) with p labelling the internal 
structure of the body. Given that the symmetric piece 
leads to a Ricci tensor, the only physically observable 
contribution will come from the quadrupolc piece Q K J — 
mpXpXp as expected. Therefore, by adding the non- 
spherical structure into its rotational part we will au- 
tomatically account for its internal quadrupole moment 
structure. Higher order multipoles do not follow this pro- 
cedure and they should be added depending on the phys- 
ical situation. 



APPENDIX D: STRESS ENERGY TENSOR FOR 
SPINNING OBJECTS 

As it was shown by Dixon a spinning particle can be 
described by the following stress energy tensor |3(| 0] , 
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14 We will show in appendix F that the spin part of the action can 
be indeed rewritten in terms of the Ricci rotation coefficients in 



a more compelling fashion. 
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The Papapetrou equations can be recovered as a con- 
sequence of Einstein equations, namely T a p.p = (I^ 15 

It can be shown also that I|D1|) is obtained from our for- 
malism. By definition the stress energy tensor is defined 
such that, 



5S 



-gT^5g aP . 



(D2) 



The variation of the action is in principle tricky due to 
the presence of the constraint e^e/„ = g^ v . Using, 



1 



0» 



S<T = -~g au g^S 9aP , 
we will therefore get, 



(D3) 
(D4) 



interesting commutation relations are those of x l a and 
spin, to leading order we have [2(j , 



Fa i X a\db 
Fa> ^a\db 



m\ 

1 



(El) 



with a = 1, 2, and db stands for Dirac brackets. 

It can be easily noticed that our expressions in I|48I5()(I 
for the potentials within the covariant condition coin- 
cides with that of 0,0,^1 U P to the acceleration piece. 
Therefore, all we have to do in order to prove the equiv- 
alence is to explicitly show that the new term derived 
from the acceleration part agrees with the extra factor 
generated by the non canonical brackets. Let us start 
with the position dynamics. It is easy to show from (|E1|) 
that the following new term, 



' f dTS IJ Sn u = ^- J d 4 xV^Tf {spm) Sg aP (D5) 

as expected. It is important to remark that, even though 
at first sight I|D1|) looks action independent, it utterly 
depends on the relationship between spin and angular- 
velocity which by itself depends on the particular La- 
grangian. 
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m-i d ( n x Si 
mi dt \ r 2 



(E2) 



will appear into the acceleration of body 1. Adding the 
piece coming from the second body it is straightforward 
to show that the extra factor is equivalently obtained by 
adding a term (in relative coordinates), 



APPENDIX E: THE EOM IN THE COVARIANT 
GAUGE: NON COMMUTATIVE ALGEBRA 

As it is known relativistic N-body EOM cannot be de- 
rived from an ordinary Lagrangian beyond the 1PN level, 
provided Lorentz invariance is preserved |42|. However, 
the latter does not follow if the position coordinates are 
not canonical variables, and this is exactly what hap- 
pens once the second class constraints in (|27|) . and the 
gauge fixing condition ip^ = 0, are imposed strongly in 
the phase space ^(J- As it is shown in [2(j the Poisson 
brackets are modified by the Dirac algebra |4fJ, |4l| . The 



J However, they do not decouple using JmJ. They can be sepa- 
rately recovered b y us ing the stress energy tensor proposed by 

Bailey and Israel plus imposing the symmetry condition 

T nv = Tv ^ by hand j 16 | 
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(E3) 



into the Lagrangian, where it is understood that wherever 
the acceleration appears in higher-order terms one sub- 
stitutes the leading order EOM. This agrees with [i"^ll3| 
and the equivalence is thus proven. In addition, it is easy 
to show there is a coordinate transformation that leads 
to a canonical algebra, to leading order [20|. 



1 



2m, 



-S a X V a . 



(E4) 



Still missing is the precession of spin. Using (|E1|) one 
can show that the EOM for spin ends up being, 

2(1 + — ) — — (n xv) x Si H j — ( n ' v i) s i- 



dt 



mi 



(E5) 



It easy to show now that the following PN shift, 

Sa^{l-\v 2 a )S a + ^V a {Va-S a ), (E6) 

which jointly with l|F4[l leads to a canonical algebra, re- 
produces the well known spin precession (see (JH2J0- The 
map in (|E4IE6|) connects the covariant choice with the 
baryonic one. The baryonic condition does not preserve 
Lorentz invariance, and an acceleration independent La- 
grangian exists 16 . As a consequence, (|51f) will not be in- 
variant under the usual linear realization of the Poincare 
group [H E2- Given that l|E4IE6|l are PN shifts it is 
also immediate to conclude that the power counting rules 
developed in this paper do not get affected by the new 
choice. 



APPENDIX F: SPIN-GRAVITON COUPLING 
REVISED 

By RPI we know that the Lagrangian must be of the 
form, 

L = -p>>u» - (Fl) 

Here we shall show that the spin part of the action can 
be rewritten as, 
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FIG. 9: A typical contribution to the one point function in 
the effective action coming from leading order mass and spin 
insertions. 



as advertised. Written this way it is clear how spin and 
gravity couple to each other, with spin playing the role 
of a 'gravitational charge' coupled to a connection of a 
spin 2 field. 



APPENDIX G: NAIVE POWER COUNTING 

Here we will schematically show the type of new in- 
sertions which are generated by divergences in the one 
point function. Let us start by calculating the effective 
action with n s spin and n rn mass leading order inser- 
tions as shown in fig. This diagram will scale as (after 
matching into NRGR for potential gravitons), 



J spin 



dT^/SuU^ 



(F2) 



with luI j the Ricci rotation coefficients. Notice this cou- 
pling is generally covariant and RPI by construction. We 
could have chosen to add spin this way into the NRGR 
Lagrangian given that S IJ can be treated as a constant 
external source (see equation H250 1. In fact, the momen- 
tum dynamics (equation (|26|l ) is also recovered following 
similar steps as we did before. It is therefore natural 
to expect that both formalisms agree which indeed fol- 
lows almost straightforwardly by definition. In terms of a 
tetrad and the Levi-Civita connection, the Ricci rotation 
coefficients can be written as [35) . 



(F3) 



Given that (|F2|) is defined on the worldline we can use 
the tetrad field transported by the particle. We will thus 
get for the spin part of the action (using l|15|)). 



= S IJ (-Tj I ^+^e a J ) = S IJ n LI , (F4) 
(It 



16 As shown by Schafer l43l . substituting the leading order EOM in 
the acceleration dependent Lagrangian of 12, 131 is also equiva- 
lent to the map to baryonic coordinates. 
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v a v m p r 
rrip" y/L r d i 



(Gl) 



where r d is introduced for dimensional reasons. By using 
NRGR power counting [l| each diagram should scale as 
/~L and therefore, 



2n„ 



1 



Lv 



2<l 



(G2) 



To consider the type of terms that can be generated 
by renormalization in the one point function, and will 
contribute to physical observables we need as a necessary 
condition d = d — 2 > 0. This however is not sufficient 
given that using Bianchi identities it can be shown 
that contraction of covariant derivatives with the full 
Riemann tensor are equivalent to derivatives of the Ricci 
tensor and henceforth f.r. removable terms. We can 
nonetheless enumerate some cases. Let us concentrate 
in logarithmic divergences first. For d = we have 
n s = n m = 1, and it is easy to see there are not any 
new terms generated by fig. The case d = 1 has 
either n s = 2, n m = or n s = 1, n m = 2. None of these 
diagrams have logarithmic divergences (for potential 
gravitons) which could generate an observable term. 
In fact, the only observable term which can be written 
down with d = 1, and either n s — 1, 2 is the first term in 
(|68p. which can be shown to be a subleading self induced 
effect. 
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For d = 2 we have either n s = 2, n m — 1 and fig. [S] which 
leads to the finite size effects we discussed in the paper, 
or n s = 1, n m = 3 which can be shown does not generate 
observable terms. For d — 3 we can have n s = 1,2,3 
plus mass insertions. Some of these diagrams will 
contribute observable terms. After matching into NRGR 
such terms start out at 0(v 10 ) for maximally rotating 
compact objects. The procedure follows systematically 
with higher order terms. 

In addition to logarithmic divergences, it is easy to 
see that diagrams like fig. |§] will also have power 
law divergences. For instance, the term O rs2 can be 
generated by diagram fig. [3] and its coefficient scale as 

I ^4 r . Assuming a cutoff of order A ~ 1/V S , we will have 

p 

C Z RS 2 ~ ~ 1/™ as expected. 

It is therefore straightforward to conclude from all we 



have seen thus far that companion induced finite size 
effect due to spin start out at 5PN (formally at 3PN) 
for maximally rotating compact objects. For the sake of 
completeness here are the terms which will contribute to 
the potential energy, 

D a D t D 1 R Pp ^S^S"'S^u\ D 2 D a R l3ppu S^S" > S^u v , 
D 2 D a R fippiJ S p °S^S v e u f} , D 2 D 1 R llvOLfi S^S° v S^u a , 
D 2 D^R^ aP S ia S afi S^u a . (G3) 

Other possible terms, like D 2 D 1 R apiUa S ia S va u p , can 
be shown to be subleading. The reasoning shown here 
can be easily extended to the case of higher order Rie- 
mann insertions with similar conclusions. 
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